Branching annihilating random walk on random regular graphs 
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The branching annihilating random walk is studied on a random graph whose sites have uniform number 
of neighbors (z). The Monte Carlo simulations in agreement with the generalized mean-field analysis indicate 
that the concentration decreases linearly with the branching rate for z > 4 while the coefficient of the linear 
term becomes zero if jz = 3. These features are described by a modified mean-field theory taking explicitly 
into consideration the probability of mutual annihilation of the parent and its offspring particles using the 
returning features of a single walker on the same graph. 
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The branching annihilating random walk (BARW) 0] 
is considered as one of the simplest models of the ex- 
tinction processes exhibiting critical behavior in differ- 
ent physical, chemical, biological and economical systems 
In these phenomena the walkers can represent do- 
main walls, vortices, defects, atoms, active sites, biolog- 
ical species or their colonies, strategies, etc. This is the 
reason why the BARWs have been extensively studied in 
the last years (for references see the papers by Cardy and 
Tauber |,|). 

In general, the walkers (henceforth particles) jump ran- 
domly on one of the neighboring sites of a lattice and 
each one can create additional particles with a branch- 
ing probability P. Furthermore, two particles annihilate 
each other if they try to share a site as consequence of the 
mentioned jump or branching events. When varying the 
branching probability a phase transition can be observed 
in the average concentration of the particles. Namely, the 
particles survive if the branching rate exceeds a critical 
value Pc, otherwise the system tends toward the absorb- 
ing state (no particles) which is independent of time. The 
transition from the active to the absorbing state belongs 
to the directed percolation (DP) universality class as well 
as the extinction processes in most of the one-component 
system 

On the lattices the BARWs are well investigated by dif- 
ferent techniques ||,|,|-|]. The Monte Carlo (MC) sim- 
ulations are extended to Sierpinski gasket by Takayasu 
and Tretyakov In the present paper we will study 
the BARW on random regular graphs characterized by 
a uniform number of joints. The joints of these graphs 
define the possible paths for the particles in the struc- 
tureless systems where the spatial position of the sites is 
no longer relevant. 

It is emphasized that the investigation of some phys- 
ical phenomena on graphs provides a more general un- 
derstanding. For example, the extension of the Mermin- 
Wagner theorem to graphs shows that the recurrence cri- 
terion for the absence of continuous symmetry break- 
ing remains valid in the graphs too In other 
words, the existence of the spontaneous magnetization 



on a graph is related to the probability of returning to 
the starting point for a single random walker on the same 
graph. The recurrence of a random walk plays also a cru- 
cial role in the BARW because a particle and its offspring 
will be annihilated when they meet. The variation of the 
distance between them can be mapped to a single walker 
problem on the same graph. As a result, if the motion 
of the parent and its offspring is not affected by other 
particles then the probability of their mutual annihila- 
tion equals to those of returning to the starting site for 
a simple random walk. 

Our investigation will be concentrated on the graphs 
consisting of N sites and each site has z joints toward 
different, randomly chosen sites (henceforth neighbors) 
excluding itself. At a given time each site can be occu- 
pied by a single particle or empty. The time evolution 
is governed by repeating the following elementary pro- 
cesses. A randomly chosen particle creates an additional 
particle on one of the neighboring site with a probabil- 
ity P or jumps to this site (with a probability 1 — P). 
In both cases, if the randomly chosen neghboring site is 
already occupied then the resident and incomer particles 
annihilate each other leaving an empty site behind. 

In the case of z = 1 the graph consists of disjoint pairs 
and the particles vanish for P > 0, while the particles 
survive on the single occupied pairs if P = 0. For z = 2 
the graph becomes a set of disjoint loops and the fea- 
ture of BARW can be described by the one-dimensional 
results [Hm^. Our analyses will be concentrated on the 
random graphs with sufficiently large N and z > 3. Lo- 
cally these graphs are similar to trees. A distance be- 
tween two sites can be introduced as the length (number 
of steps) of the shortest path joining them. The average 
distance between two sites increases logarithmically with 
the number of sites for large N [jl^ . 

The stationary state of this system is characterized by 
the average concentration of walkers (c) that will be de- 
termined by using different methods. For a locally tree 
like structure the generalization of the one-dimensional 
dynamical cluster technique is straightforward ]l3[ |. In 
this case the particle distribution is described by the 
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configuration probabilities pk{ni, ■ ■ ■ ,nk) {rii = or 1) 
on the clusters of the neighboring k sites. Here we as- 
sume that these quantities satisfy some symmetry (trans- 
lation, rotation, reflection) and compatibility relations. 
The one-point configuration probabilities are directly re- 
lated to the average concentrations, namely, pi{l) = c 
and pi{0) — I — c. Introducing an additional parameter 
q, the two-point configuration probabilities are given as 
P2{1, 1) = q, P2(l, 0) = p2(0, 1) = c - g and p2(0, 0) = 
1 — 2c + q. Further parameters are required for k > 2. 

In the present case the time variation of pk can be 
expressed by the terms of pk and Pk+i- For example. 



Pl(l) = (1 - P)pi(l) - P2{1, 1) +P2{0, 1) 



(1) 



where we have summed the contribution of all the el- 
ementary processes mentioned above. Notice that this 
equation is satisfied by the absorbing state (c — 0). 
At the level of one-point approximation we assume that 
P2{ni,n2) — Pi{ni)pi{n2) ■ In this case the non-trivial 
stationary solution of Eq. (|]) obeys a simple form. 



(2) 



independent of z. At the level of fc-point approximation 
the corresponding set of equations is solved by using the 
Bayesian relations (pfc+is are approximated by the prod- 
uct oipk terms) |13 . In the two-point approximation the 
straightforward calculation gives the following stationary 
solution: 



^(2rt^p.2(^-2)-(z-3)P 



4(z - 1) - 2zP 2P2 



(3) 



This result refers to the absence of pair correlations in 
the limit z — > cxd as well as at P = 1 for any values of 
z. At higher levels the stationary solutions are evaluated 
numerically. 

In order to check these results we have performed MC 
simulations on random graphs with N = 500, 000 sites 
varying the branching rate P for z — 4: and 3. The simu- 
lations are started from a randomly half-filled graph and 
the concentration is obtained by averaging over 10'* MC 
steps per particles after some thermalization. 

Figure |l| compares the MC data to the prediction of 
fc-point approximations for z = 4. Here the results of 
3- and 4-point approximations are omitted because their 
deviation from c'^^^-' is comparable to the line thickness. 
In this case c = AP when P ^ 0. The A coefficients 
obtained by MC simulation and 5-point approximation 
are slightly different, namely, A^*^"^) = 0.250(2) and 
^(5p) ^ 0.278(1). 

In the light of pair approximation [see Eq. (||)] better 
and better agreement is expected when increasing z. For 
z — 3, however, significant differences can be observed 
between the MC results and the prediction of fc-point 
approximations when P — > 0. To magnify the discrep- 
ancy between the two methods the P-dependence of the 



average concentration is illustrated in a log-log plot. As 
shown in Fig. || the MC data refer to a quadratic behav- 
ior for small P values, meanwhile the fc-point approxima- 
tions predict linear behavior with a coefficient decreasing 
when k is increased. 




FIG. 1. The average concentration of particles as a func- 
tion of branching rate for z = 4. The symbols represent MC 
data, the solid curves comes from the 1-, 2- and 5-point ap- 
proximations (from top to bottom) on clusters indicated at 
the top. 




FIG. 2. Log-log plot of the average concentration of parti- 
cles vs P if z — 3. The solid curves indicate the prediction 
of generalized mean-field methods at the levels of 1-, 2-, 4-, 
and 6-point approximations (from top to bottom). The fitted 
curve (dashed line) on the MC data (open diamonds) shows 
quadratic behavior. 

Here it is worth mentioning that the MC data have 
remained unchanged within the statistical error (compa- 
rable to the line thickness) when the random graph was 
generated in a different way. This investigation was mo- 
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tivated by the small word model suggested by Watts and 
Strogatz In this case the points with the first and 
second joints form a single loop and the third ones are 
chosen randomly. 

The trend in the prediction of fc-point approximations 
implies the possibility that the coefficient of the linear 
term goes to zero in the limit k oo. Now wc introduce 
a modified mean-field theory taking explicitly into con- 
sideration the mutual annihilation of the parent and his 
offspring as mentioned above. It will be shown that the 
vanishment of the linear term is directly related to fact 
that half of the branching process yields mutual annihi- 
lation. 

In this description the R{z) probability of this mutual 
annihilation process is approximated with the probability 
of returning to the starting site for a single random walker 
on the same graph. This simplification is exact in the 
limits P and c — > when the random walks of the parent 
and his offspring are practically not affected by other 
particles. The appearance of additional offsprings during 
the recurrence is considered a second order process whose 
effect will be discussed below. 

The main advance of this approach is that we can use 
the exact results obtained for the simple random walk on 
the Bethe lattices which are locally similar to the present 
random graphs. Hughes and Sahimi [|5| and Cassi 
have obtained that R{z) = l/(z — 1) and the average 
number of steps to return to the starting site is r(z) = 
2{z — l)/{z — 2). The value of t{z) indicates that the 
particle visits only a few sites before its recurrence and 
during this short period it is not capable to distinguish 
the Bethe lattices from the present random graphs due 
to the absence (low probability) of loops. In fact, this is 
the reason why our analysis is restricted to large N. 

As a result of the mutual annihilation, R{z) portion of 
the branching process can be considered as spontaneous 
annihilations. On the other hand, 1 — R{z) portion of 
the branching events results in "independent" particles. 
These features can be easily built into a modified mean- 
field theory which obeys the following form: 

= - (1 - P)pi(l) -P2(l, 1) - Pi?(z)p2(l,0) 

+ [l-PR{z)]p2{0,l). (4) 

The two-point configuration probabilities are also af- 
fected by the random walk itself, because a given particle 
leaves an empty site behind when it steps to one of the 
neighboring sites. Consequently, the value of p2(l,0) is 
larger than those predicted by the mean-field theory. The 
above mentioned techniques confirm that this correlation 
can be well described by a simple parameter defined as 



P2(1,1) = QWpi(1)pi(1) 



(5) 



where Q{z) < 1 and the remaining two-point configu- 
ration probabilities are determined by the compatibility 
conditions. The value of Q{z) is related to the asymp- 
totic time-dependence for the annihilating random walk 



{P = 0) when the particle concentration c{t) decreases 
monotonously. From Eqs. (Q) and ^ we obtain that 

when t ^ oo. Figure || compares this result with the MC 
data as well as with the prediction of fc-point approxi- 
mations obtained by numerical integration of the corre- 
sponding set of differential equations for z = 3 and 4. 
In the MC simulations the system is started from a half- 
filled random initial state for N — 10^ and the data are 
averaged over 50 runs. The classical mean- field (1-point) 
approximation corresponds to the choice oi Q{z) — 1. 
From the fc-point approximations, however, we can de- 
duce smaller Q{z) values. For example, in the pair ap- 
proximation Q(2P)(4) = 3/4 and Q'^^p^S) = 2/3. These 
estimations can be improved if we choose larger and 
larger clusters as illustrated in Fig. We have obtained 
Q(5p)(4) = 0.692(1) for z = 4 and Q^p\3) = 0.534(1) for 
z — 3 using 5- and 6-point approximations respectively. 
At the same time, the numerical fitting to the MC data 
results in Q(^'=')(4) = 0.67(2) and Q^^^Hi) = 0.50(2). 
Notice that the predictions of fc-point approximations 
tend slowly toward the MC values. According to the nu- 
merical integration in fc'-point approximations the asymp- 
totic behavior is independent of the initial concentration 
c(0). 




10" 
f [MCS] 

FIG. 3. Time-dependence of average concentration multi- 
plied by time for P = and c(0) = 1/2. The open squares 
(diamonds) represent MC data for 2 = 3 (z = 4). The 
upper (lower) solid curve indicates the prediction of 6-point 
(5-point) approximation for z = 3 (z = 4). The dashed line 
shows the prediction of classical mean-field approximation. 
The fitted asymptotic functions are represented by dotted 
lines. 

Using the expression (H) the nontrivial stationary so- 
lution of Eq. (^) obeys the following simple form: 



c = P 



1 - 2R{z) 



2[1 - PR{z)]Q{z) 



(7) 
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Evidently, beside this expression Eq. has a trivial 
solution (c = 0) which remains the only one for R > 1/2. 
In the limit P ^ the leading term of Eq. can be 
expressed as c = P[l — 2R{z)]/2Q{z). This prediction 
coincides with the above mentioned MC result for z = 4 
if Q(*^^)(4) is substituted for Q{z). 

According to Eq. the average concentration be- 
comes zero for z = 3. This is a consequence of the fact 
that the rare branching processes do not modify the av- 
erage number of walkers because these events result in 
zero or two walkers with the same probability on a longer 
time scale. There exists, however, a second order term 
neglected above which is responsible for the quadratic be- 
havior. Namely, during the recurrence one of the (parent 
and offspring) particles can create an additional walker 
with a probability proportional to P'^R(z)[t{z) — 1] and 
this event reduces the rate of spontaneous annihilation 
[PR{z)] in Eqs. (^) and (0). Unfortunately, this rough 
approach can not reproduce accurately the coefficient of 
the quadratic term because of the simplification we used 
in the derivation of Eq. (^. 

The present analysis indicates that the probability of 
recurrence of a single walker, the long time behavior of 
the time-dependent concentration for the annihilating 
random walks and the stationary concentration for the 
BARW in the limit P ^ are strongly related to each 
other. The modified mean-field theory suggests a way to 
describe this relation on random regular graphs for suffi- 
ciently large number of sites. The relation is based on two 
simple conditions. Namely, R{z) < 1/2 and t{z) is finite. 
These conditions are also fulfilled on the cubic lattices for 
dimensions d > 4 [T^ where the rigorous analysis predicts 
mean- field type behavior Q]. Conversely, the modified 
mean-field theory is not adequate for lower dimensions 
{d < 3) because r = oo. It is expected, however, that the 
present analysis can be extended to the investigation of 
BARW on such random graphs characterized with differ- 
ent probability distributions of connectivity |l^Jl^] . The 
present random graphs for z = 3 represent a particular 
situation where the coefficient of the linear term vanishes 
as a consequence of i?(3) = 1/2. The solutions of the 
modified mean-field theory imply the the possibility of 
another interesting situation for those above mentioned 
random graphs [^14 l^Jl^ where the average value of 
R is larger than 1/2. In this case the particles can sur- 
vive if the branching rate exceeds a treshold value as it 
happens on the one- and two-dimensional lattices. 
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